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ABSTRACT 



Context. Detailed hydrodynamic modeling explained several features of a fragment of the core of a Coronal Mass Ejection observed with 
SoHO/UVCS at 1.7 Rq on 12 December 1997, but some questions remained unsolved. 

Aims. We investigate the role of the magnetic fields in the thermal insulation and the expansion of an ejected fragment (cloud) traveling 
upwards in the outer corona. 

Methods. We perform MHD simulations including the effects of thermal conduction and radiative losses of a dense spherical or cylindrical 
cloud launched upwards in the outer corona, with various assumptions on the strength and topology of the ambient magnetic field; we also 
consider the case of a cylindrical cloud with an internal magnetic field component along its axis. 

Results. We find that a weak ambient magnetic field (J3 ~ 20) with open topology provides both significant thermal insulation and large 
expansion. The cylindrical cloud expands more than the spherical one. 
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1. Introduction 

The study of Coronal Mass Ejections (CME) has received 
a strong impulse since the la unch of the Solar Heliosp heric 



a strong impulse since the launch or the solar Heliospheric 
Observatory (SoHO) in 1995 dDomingo & Polandll 19881) . The 



Ultraviolet Coro nal Spectrometer (UVCS) on board SoHO 
( Kohl et al. 1995b monitored the evolution of many CMEs and 
brought detailed spectroscopic information on plasma in CMEs 
during their propagation in the outer corona. 

Ciaravella et al. d2000h and ICiaravella et al.l 
lyzed in detail a CME occurred on 12 December 



(1200 lb ana- 
1997. 

Ciaravella et al.ld2001l) addressed more specifically the features 
observed at a distance of about 2 Rq by SoHO/UVCS, and hy- 
drodynamic (HD) simulations of a cloud expelled out of the 
corona were performed to investigate which features of a frag- 
ment of the CME core observed with SoHO/UVCS could be 
explained in terms of purely HD effects. That model assumed 
the CME fragment to be already formed and launched with an 
initial impulse, and studied its evolution while traveling from 
the height of the EIT observation (~ 1.2 Rq) to that of the 
UVCS observation (~ 1.7 Rq). The HD model was able to 
reproduce various features observed in SoHO/UVCS UV lines 
(C III and O VI). The fragment was constrained to be initially 
denser than its surrounding with a density contrast dn/no ~ 3 
and the background density to be ~ 10 s cm~ 3 . 



From the results of HD modeling. ICiaravella et alj (2001) 
argued that the thermal insulation of the core fragment is re- 
quired to explain its evolution and to fit the UVCS observa- 
tion. However, two main questions remained open. First, the 
thermal insulation needs to be physically motivated, since the 
thermal conduction is very efficient in the corona; in particu- 
lar numerical simulations showed that the CME fragment ther- 
malizes with the hot atmosphere on a very short time scale un- 
less the thermal conduction is somehow suppressed. Second, 
the CME is observed to expand by a factor 3 to 4 in diame- 
ter (ICiaravella et al . 2000) traveling from the low corona to the 
height of the UVCS observation, while the hydrodynamic sim- 
ulations lead to a smaller expansion of the fragment. Although 
the modeling suggests that very elongated structures may ex- 
pand more, the amount of expansion needs further investiga- 
tion. 

Here, we extend the modeling study described in Ciaravella 
et al. (2001) and investigate whether magnetic fields may have 
a crucial role to answer the questions left open in the previ- 
ous study, i.e. the thermal insulation and/or the amount of the 
expansion of the moving CME fragment. The magnetic field 
is able to funnel the thermal conduction, and it may lead to 
a strong thermal insulation of the fragment depending on the 
magnetic field topology. F or this purpose, we revi sit the model- 
ing approach described in ICiaravella et al.l d200lb now includ- 
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ing an ambient magnetic field and using a full MHD modeling 
that includes non-isotropic thermal conduction, to understand 
in which conditions the expelled fragment can be thermally in- 
sulated or may expand more. We study the evolution of the 
moving fragment for different topologies and intensities of the 
ambient coronal magnetic field, and eventually consider also a 
fragment with an internal magnetic field, therefore accounting 
for the ejection of magnetic flux. In section |2] we describe the 
model, section [3] presents the simulations and in section [4] we 
discuss the results. 

2. The model 

Our scope is to study the evolution of an ejected CME core 
fragment as it travels through the high corona, extending the 
modeling described in Section 5.2 of Ciaravella et al. (2001). 
Again, we consider a simplified model of a plasma cloud, 
denser than the surrounding atmosphere and moving upwards 
in a stratified atmosphere with parameters that lead to the 
conditions observed at 1.7 Rq with UVCS. At variance with 



Ciaravella et al.l (120011) now we consider the presence of an am- 
bient magnetic field. W e model the fragment while it travels in 
a [3 i5 1 environment ( Garvl 200ll) . where we cannot neglect 
the magnetic pressure, nor the thermal one, and we include the 
effect of the solar gravity, the radiative losses from an optically 
thin plasma, a coronal heating (needed to keep the unperturbed 
atmosphere in hydrostatic conditions), and the anisotropic ther- 
mal conduction in the presence of a magnetic field. 

For continuity with the work of Ciaravella et al. (2001), we 
will first report on the evolution of a plasma cloud moving in 
various configurations of ambient magnetic field and eventu- 
ally consider the case of a cloud of magnetic flux and plasma, 
i.e. a "magnetized" cloud, with an internal magnetic field. 

We solve the ideal MHD equations (here in CGS conserva- 
tive form): 



| + V.( P .) = 



(1) 



dpv „ . , „ (V x B) x B 

+ V • (pvv) = -Vp + ^ + pg (2) 



du 
dt 



V • [(it + p)v] =pgv- n 2 P(T) + H - V ■ F c (3) 



dB 

-— = V x (v x B) 

dt v ; 

1 2 
u = —pv + E 

P = h-i)E 



(4) 

(5) 
(6) 



where t is the time, p is the mass density, n the number density, 
p the thermal pressure, T the temperature, v the plasma flow 
speed, E the internal energy, u the total energy (internal plus 
kinetic), 7 the adiabatic index, g the gravity acceleration, B 
the magn etic field, P(T) the radia tive losses per unit emission 
measure (Raymond & Smith 19771) . H a space-dependent heat- 
ing function equal to the radiative losses in the initial conditions 



and F c is the conductive flux according to Spitzer! dl962). We 
neglect resistivity effects. 

We solve numerically the equations with the MHD version 
of the advanced parallel FLASH code, basically developed by 
the ASC / Alliance Center for Astrop hysical Thermonuclear 
Flashes in Chicago jFrvxell et al.ll2000[). with Ad aptive Mesh 
Refinement PARAMESH dMacNeice et alJl2000h . The MHD 
equations are solved using the numerical scheme HLLE 



(flatten Lax Van Leer Einfeldt) proposed by Einfeldt (1988). 
We implemented a FLASH mod ule to model t he anisotropic 
thermal conduction according to Spitzer (1962). The thermal 
flux is locally split into two components, along and across the 
magnetic field, and is given by: 



F c = -/s,|(VT)|| -kl(VT)_l 



(7) 



where (VT)m and (VT)j_ are the thermal gradients along and 
across the magnetic field, and Kit and k± the thermal conduc- 
tion coefficients alon g and across the magnetic field lines, re- 
spectively, given by Spitzer ( 1962): 



S T e20 



3/2 



{k b Tf' 2 k b 



1/2, 



'Zln(A) 



8{irm l k b ) 1 / 2 n 2 Z 2 e 2 c 2 ln(A) 
ZB 2 T X I 2 



(8) 



(9) 



where c is the light speed, k b is the Boltzmann constant, m e 
and rrii are respectively the electron mass and the average ion 
mass, e is the electron charge, Z is the average atomic number, 
8t and e parameters depends on the chemical composition, and 
in a proton-electron plasma: 8t — 0.225 and e = 0.419; ln(A) 
is the Coulomb logarithm which is (for the coronal plasma) 
Ln(A) ~ 20. 

The unperturbed corona is assumed to be in magnetohydro- 
static conditions, isothermal (T = 1.5 x 10 6 K) and stratified 



by gravity with a number density uq = 10 



at the initial 



height of the cloud centre. The cloud, assumed in pressure equi- 
librium with the atmosphere, is denser (jicme = ^iatm) and 
colder than its surrounding and its center is a height of ~ 0.15 
Rq above the photosphere . It is also stratified to p reserve initial 
isobaric equilibrium. As in lCiaravella et al. (2001), the cloud at 
t = s is assumed circular, which in a 2-D cartesian geometry 
is equivalent to a cylindrical cloud with infinite y-extension in 
a proper 3-D domain. 

For all the simulations, the radius of the cloud is 5 x 10 9 cm 
and the initial upward velocity of the cloud is 400 km/s, super- 
sonic but smaller than the escape velocity (~ 600 km/s), and 
large enough to reach the height of the UVCS observation at 
1.7 Rq. The simulations describe the evolution of the plasma 
for a time interval of 3000 s, i.e. the time taken by the cloud 
to reach the UVCS observation height. The coordinate system 
is 2-D cartesian (x,z), the vertical direction is the z coordinate. 
The computational domain extends from the solar photosphere 
(z = cm) to « Rq (z = 7 x 10 10 cm) in the z direction and to 
f» 0.4i? Q in the x direction (x = cm to x = 3 x 10 10 cm). The 
symmetry of the system along the x = cm axis allows us to 
model half spatial domain. The adaptive mesh algorithm yields 
an effective resolution of ~ 3 x 10 8 cm (16 zones per cloud 
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radius). Reflection boundary conditions are imposed along the 
x = cm axis consistently with the symmetry. At the upper 
and external boundaries we set zero-gradient (outflow) bound- 
ary conditions. Fixed values are imposed at the lower boundary 
(z = cm). We use 2-D geometry for most of our modeling. 

From this setup we perform a basic set of simulations with 
different ambient coronal magnetic fields, B, in which the 
cloud evolves, i.e. (a) with no magnetic field, (b) with a dipolar 
magnetic field on the equatorial plane {closed field), (c) with a 
dipolar magnetic field on the polar axis (open field). Since we 
follow the cloud evolution much later than the launch and far 
from the origin site, our ambi ent magneti c field is relatively 
weak, as of the outer corona dGarvl 120011) . In particular, we 
assume that the dipole lays in the center of the Sun and its 
strength is tuned so to have j3 ~ 25 at the cloud height, as basic 
value. Figure[T]shows the initial conditions for the simulations 
with an ambient dipolar field {closed and open). In such con- 
figurations, the highest (3 rs 40 is at the base of the atmosphere 
(where the pressure is higher), (3 varies from 30 to 20 inside the 
initial cloud, and then it settles to (3 ~ 15 from a height of 1.5 
i?© up. Simulations with more strong magnetic fields have also 
been considered (see below), and there the fractional variation 
of (3 is the same throughout the spatial domain. 

Other simulations are also performed: we repeat simula- 
tions with no magnetic field and with open dipole magnetic 
field for a sph erical cloud (ins t ead of cylindrical), the same 
as presented in ICiaravella et al.1 (120011) : we make some simu- 
lations with a stronger ambient magnetic field (i.e. (3 ~ 0.25 
at the height of the centre of the initial cloud), a magnetic- 
field-dominated initial regime. To account for the ejection of 
both magnetic flux and plasma we replicate the simulation with 
the open dipole ambient magnetic field for an initially magne- 
tized cylindrical cloud, i.e. an additional magnetic field compo- 
nent inside the cloud along y (i.e. along the axis of the cloud), 
as shown in Fig. [2] We consider two different values of B y : 
B y = 0.2 G and B y = 1 G. With B y = 0.2 G, the cloud 
parameters are left unchanged, therefore creating a small over- 
pressure on the ambient medium. With B y = 1 G, the cloud 
temperature is set equal to 40,000 K (i.e. much lower than the 
reference value) to preserve the initial condition of isobaric 
cloud. 

Table[T]lists the relevant simulations presented in this work, 
identified by the presence or the absence of magnetic field, the 
magnetic field topology, the value of (3 at the initial height of 
the cloud centre, the cloud shape, the magnetic field inside the 
cloud, the presence or absence of thermal conduction, the ge- 
ometry of the computational domain. 

The simulation with open magnetic field and a spherical 
cloud requires 3-D modeling and we extend the cartesian do- 
main also in the y-direction to y = 3 x 10 10 cm. We have per- 
formed these simulations on a Cluster Linux EXADRON with 
24 Opteron 250 AMD processors and on a IBM SP Cluster 
with 512 IBM Power5 processors. The 2-D simulations have re- 
quired ~ 300 hours of computational time and the 3-D ~ 3000 
hours. 
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Fig. 1. Magnetic field lines and color maps of the density con- 
trast (p/po) for two different initial configurations of mag- 
netic field: closed field (upper panel, those starting with MC 
in TableHJ and open field (lower panel, those starting with MO 
in Table [T). 



3. Results 

3.1. Cylindrical cloud 

3.1 .1 . Simulations with no magnetic field 

We start the simulation with a cylindrical cloud in the envi- 
ronment without magnetic field as the basic model (case HC, 
abbreviation f or Hydrodynamic mode l, Cylindrical cloud, in 
Table [TJ. As in Ciaravella et al. ( 200 lb the thermal conduction 
is set to zero in this simulation. Figure [3] shows color maps of 
the density contrast (ratio of the density to the density of the 
static atmosphere p/po) and the temperature at times t = 500 
s, 1500 s,and 3000 s, i.e. at two intermediate times and at the 
final time. 

The cloud initial velocity is higher than the local sound 
speed (eg ~ 200 km/s) and a shock front soon departs radi- 
ally from the cloud. While moving upwards, the cloud expands 
and dynamic instabilities develop at its boundary, changing its 
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Table 1. Numerical simulations 



Name 


Magn. 


B dipole 






Cloud 


Cloud 


Thermal 


Geometry 




field (B) 


topology 






shape 


internal B v 


conduction 




HC 


No 








Cylinder 




No 


2-D cartesian (x, z) 


HS 


No 








Sphere 




No 


2-D cylindrical (r, z) 


MCLC 


Yes 


Closed 


Low (~ 


0.25) 


Cylinder 




Yes 


2-D cartesian (x, z) 


MOLC 


Yes 


Open 


Low (~ 


0.25) 


Cylinder 




Yes 


2-D cartesian (x, z) 


MCHC 


Yes 


Closed 


High (r 


- 25) 


Cylinder 




Yes 


2-D cartesian (x, z) 


MOHC 


Yes 


Open 


High (r 


- 25) 


Cylinder 




Yes 


2-D cartesian (x, z) 


MOHS 


Yes 


Open 


High (r 


- 25) 


Sphere 




Yes 


3-D cartesian (x, y, z) 


MOHCL 


Yes 


Open 


High (r 


- 25) 


Cylinder 


Low (0.2 G) 


Yes 


2-D cartesian (x, z) 


MOHCH 


Yes 


Open 


High (r 


- 25) 


Cylinder 


High (1 G) 


Yes 


2-D cartesian (x, z) 
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moves upwards, we define an expansion factor (F e ) as the ratio 
between the maximum width of the cloud (in the ^-direction) 
and its initial radius, which can be c ompared to the factor 3-4 
estimated bv lCiaravella et al" ( 2000l) . At the end the expansion 
stops and the pressure equilibrium between the cloud and its 
surroundings is recovered. We obtain F e « 2.9 at t=3000 s for 
this simulation. 

We do not prese nt here hydrodynamic s imulations with the 
thermal conduction. Iciaravella et all ( 2001 ) showed that if the 
thermal conduction were not suppressed, the core would mostly 
evaporate because of heating by the hot surrounding corona and 
would shrink to a very small cold knot. 



Fig. 2. Sketch of the initial configuration of the simulations 
with an open dipolar ambient magnetic field and an initially 
magnetized cloud (MOHCL and MOHCH in TableQ). The am- 
bient dipolar magnetic field is shown in the x-z plane (see lower 
panel of Fig.Q). The additional magnetic field component is set 
inside the cloud (the long cylinder) along the y-direction. 

shape and forming small scale structures, departing from it. At 
the final time the cloud has evolved into a cold extended core 
with two thin tails and has mostly lost memory of its initial 
circular shape. The core becomes colder and denser during the 
evolution because of the radiative losses and of the absence of 
thermal conduction. Since the cloud is stratified, thermal in- 
stabilities first occur in the lower (and denser) arc-shaped part, 
and at the end they involve the whole cloud. Because of the 
relative motion between the cloud and the atmosphere, Kelvin- 
Helm holtz and Rayleigh-Taylo r instabilities develop in a time- 
scale (Chen & Lvkoudislll972h : 



trhrt 



_J_cs_ 

7 - 1 9 



10 3 . 



(10) 



where c s is the sound speed. The motion of the cloud is mostly 
ballistic, although it loses part of its mechanical energy (i.e. 
kinetic plus potential gravitational) in the hydrodynamic inter- 
action with the atmosphere, and at 3000 s the cloud has practi- 
cally stopped. To quantify the expansion of the cloud while it 



3.1 .2. Simulations with magnetic field 

We now report on results of simulations with an initially cylin- 
drical cloud, with non-zero ambient magnetic field, and with 
the thermal conduction strongly effective along the magnetic 
field lines (MCLC, abbreviation for Magnetohydrodynamic 
model, Closed field, Low (3, Cylindrical cloud, MOLC, MCHC, 
and MOHC in TableQ). 

In the low j3 simulations (MCLC, MOLC) the cloud does 
not show significant evolution, because the magnetic field is 
too strong to be perturbed by the moving cloud. The strong 
magnetic tension does not allow the cloud to move upwards, 
either. As an example, figure |4] shows the case of the closed 
dipolar field (MCLC) at t = 500 s, with (3 ~ 0.25 at the initial 
height of the cloud centre. 

In the high (3 ambient medium, instead, the cloud moves 
upwards and expands in the outer corona. Figure [5] shows the 
evolution of the density contrast and the temperature with color 
maps at t = 500 s, 1500 s and 3000 s for the simulation with 
the closed dipolar field for (3 ~ 25 at the initial cloud height. 
Since the initial cloud velocity (400 km/s) is higher than both 
the sound speed and the Alfven speed (va ~ 10 7 cm/s), a fast 
MHD shock propagates radially from the cloud. During the 
evolution no hydrodynamic instability develops because of the 
suppression by the thermal conduction and the magnetic field. 
In fact, they are effective on time scales smaller than tkhrt 
(Eq.fTOli, i.e. of the order of, respectively: 



T c 



21nk b L 2 
2fc T 2 - 5 



100 s 



(11) 
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Fig. 3. Color maps of temperature (left of each panel) and den- 
sity contrast, p/po (right of each panel) at times t = 500 s, 
1500 s and 3000 s, for the simulation with no magnetic field 
(and no thermal conduction, HC in Table[TJ. 



— - 100 s 

VA 



(12) 



where fc T 2 5 = k\\ and for typical coronal density with k = 
9.2 x 10~ 7 (c.g.s. units, from EqJH]), va the Alfven speed, L a 
characteristic length scale, which we have conservatively taken 
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Fig. 4. Magnetic field lines and color maps of temperature and 
density contrast at time t = 500 s for the simulation with a 
strong closed magnetic field (MCLC in Table [T). The magnetic 
field strength is w 2.7 G at z ~ 1.5 x 10 10 cm along the x=0 
axis (high density field lines) and no field line was drawn in the 
region where it becomes less than ~ 1 G (z > 5 x 10 10 cm). 



as half of the cloud radius. The cloud drags the magnetic field 
lines which gradually envelope the cloud and inhibit thermal 
exchanges between the cloud and the surrounding medium. The 
consequent strong thermal insulation leaves the plasma free to 
cool down to temperatures below T rs 10 5 6 K at which it emits 
the UV lines of interest. The upward motion of the cloud is 
braked by the restoring forces of the magnetic field; only the 
upper part of the cloud reaches the UVCS observation height. 
The cloud expansion is limited by the confining effect of the 
magnetic field, which becomes up to 5 times stronger just out- 
side the cloud (J3 ~ 0.1) than inside it. At the end, the expan- 
sion factor of the cloud is F e « 2.1 at t ~ 3000 s (i.e. ~ 70% 
of the reference case). 

The cloud travels much longer distances when it moves in 
an open ambient magnetic field. Figure [6] shows the evolution 
of the density contrast and of the temperature at t = 500 s, 1500 
s and 3000 s for the related simulation. Again, a fast MHD 
shock propagates upward radially from the cloud and no hy- 
drodynamic instability develops. At variance with the closed 
field case, in the open field the cloud moves ballistically to 
the outer corona. The cloud is thermally insulated laterally, 
because the initial field direction is mostly maintained during 
the evolution. In the upper corona, the magnetic field is weak 
and, therefore, perturbed by the rising cloud (at t > 2000 s). A 
strong magnetic field horizontal component develops ahead of 
the cloud and inhibits thermal conduction with the region above 
the cloud. As the magnetic field in the cloud is frozen and the 
cloud moves fast with respect to the Alfven scale times, it drags 
the magnetic field lines, and inversely-directed downward mag- 
netic field components and current sheets are produced near the 
cloud. The open magnetic field naturally favours the cloud ex- 
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Fig. 5. Magnetic field lines and color maps of temperature and Fig. 6. Magnetic field lines and color maps of temperature and 
density contrast at t=500 s, 1500 s and 3000 s for the Simula- density contrast at t=500 s, 1500 s and 3000 s for the simula- 
tion with weak closed magnetic field (MCHC in Table [TJ. For tion with weak open magnetic field (MOHC in Table [TJ. The 
reference, at t=500 s the strength of the magnetic field ranges strength of the magnetic field ranges from % 0.03 G to « 0.5 
between % 0.05 G (in the upper part of the domain where no G. 
magnetic field lines are drawn) and « 0.7 G (just above the 
cloud where the lines are very dense). 

3.2. Spherical cloud 

pansion and, in fact, the final expansion factor is F e « 3.9 at We now present results for the case of a spherical cloud. We 

t ~ 3000 s (i.e. ~ 35% more than the reference case). compare the null-field evolution (the one presented with the 
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highest detail in ICiaravella et al.1 d200ll) . HS in Table Q) with 
the evolution of the cloud in the weak open dipole field (j3 ~ 25 
at the initial cloud height, MOHS in Table [TJ, performed with 
a fully 3-D simulation. Figure [7] shows the density contrast for 
both the simulations. 
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Fig. 7. Magnetic field lines and color maps (on the section 
across the cloud axis) of density contrast at t=500 s, 1500 s and 
3000 s for the simulations of a spherical cloud without mag- 
netic field (left of each panel, HS in Table [TJ and with a weak 
open magnetic field (right of each panel, MOHS in Table [TJ. 
The strength of the magnetic field ranges from w 0.03 G to 
*0.5G. 

In the HD simulation (HS), the cloud evolves into a thin 
dense shell-like structure, with very irregular boundaries, due 
to the hydrodynamic instabilities, that can develop in the ab- 



sence of thermal conduction. At the end of the evolution, the 
expansion factor of the cloud is F e as 2.7 at t ~ 3000 s. Since 
the volume expansion of the spherical cloud scales as r 3 , the 
corresponding projected expansion is less significant than the 
projected expansion of the cylindrical cloud, whose volume 
expansion scales as r 2 . In the MHD simulation (MOHS), the 
more effective expansion makes the cloud apparently less con- 
spicuous than in the case of the cylindrical cloud (Figj6]). The 
expansion factor is F e = 2.8 at t ~ 3000 s, ~ 4% more than 
the HD case, and ~ 70% of the expansion of the cylindrical 
cloud (MOHC, see Section[3X2l>. 

3.3. Magnetized cloud 

As mentioned in Section 12 we consider two simulations 
of magnetized clouds, one with a weak (0.2 G) inter- 
nal £? y -component (i.e. along the axis of the cylindrical 
cloud, as sketched in Fig. 12 MOHCL (abbreviation for 
Magnetohydrodynamic model, Open field, High f3, Cylindrical 
cloud, Low internal magnetic field) in Table [TJ the other with 
a strong (1 G) S^-component (MOHCH) in Table [TJ We will 
not comment much on the latter case, because we found that the 
whole cloud remains very cool all over the computed evolution, 
much cooler than the temperature of formation of the relevant 
UVCS lines. This case therefore seems to be unrealistic. 

Figure[8]shows magnetic field, density contrast and temper- 
ature distributions over the cloud cross-section for the simula- 
tion with the weaker magnetic field component inside the cloud 
(MOHCL). The global evolution does not change significantly 
with respect to the case of no internal magnetic field (MOHC, 
Section [3. 1.2b . and the morphology of the cloud evolves simi- 
larly to that of Figure [6] Therefore, the presence of an internal 
magnetic field component does not influence much the evolu- 
tion of the cloud, in spite of the initial cloud overpressure. The 
latter is only a small cloud rapidly absorbed by the corona in 
t ~ L/ c§ ~ 250s, where L = 5 x 10 9 cm is the cloud size. 
The contours on the right side of the panels of Fig.[8]show that 
the By component initially set up inside the cloud becomes 
weaker and weaker as the cloud moves upwards and expands. 

The cloud thermal insulation is very good in this simula- 
tion, as indicated by the growing cold core, even more extended 
than in the other simulations. The magnetic field component 
along the axis of a cylinder infinitely extending horizontally 
further inhibit thermal exchanges with the surrounding corona. 
In a more realistic configuration of an elongated cloud with 
an internal magnetic field which bends downwards and even- 
tually connects to the photosphere (e.g. a flux rope), some heat 
would be conducted to the footpoints, similar to the 1-D con- 
duction case described in Ciaravella et al. (2001). We estimate 
that the cloud would thermalize in r c ~ 1000 s (see EqJTTb. 
if the lenght scale of the cloud is the same as its height, say 
h ~ 10 10 cm, above the photosphere. The cloud expansion fac- 
tor for this case of magnetized cloud is F e w 3.6 at t ~ 3000 s. 

4. Discussion and conclusions 

This work is devoted to studying the role of the magnetic fields 
in the evolution of a fragment of a CME core traveling upwards 
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expla ined with a purely hydrodynamic model dCiaravella et al 



2001). In the simulations presented here, we basically consider 



Fig. 8. Magnetic field lines (x-z plane, left of each panel), color 
maps of temperature {left) and density contrast {right), and con- 
tours of the magnetic field y-component {right) at time t=500 s, 
1500 s and 3000 s for the simulation with a weak open ambient 
field and a weak magnetic field y-component inside the cloud 
(MOHCL in Table [Q). The contours of B y are spaced by 0.02 
G to a maximum of 0.2 G. The strength of the magnetic field 
on the left of the panels ranges from w 0.03 G to w 0.5 G. 



in the high solar corona. We address the late stage evolution of 
the cloud, in which it travels in a weakly magnetized atmo- 
sphere {(3 <; 1). Our simulations show that the evolution in 
the late stage could be described better than in the early stage. 
In particular, we focus on the thermal insulation of the cloud 
and on its degree of expansion; both aspects were not well 



infinitely long cylindrical clouds, and we first investigate which 
configuration (i.e. strength and topology) of the ambient mag- 
netic field could favor better the thermal insulation and the ex- 
pansion of the cloud simultaneously. 

Before investigating the cloud enclosed by magnetic fields, 
we screen out those ambient magnetic field configurations 
which strongly brake the upward motion of the cloud, i.e. pro- 
vide magnetic confinement. Since the plasma is frozen to the 
field, the cloud moving upward drags the ambient magnetic 
field. The strong magnetic tension then acts against the cloud 
expansion and motion. Our simulations show that the cloud is 
strongly braked in an atmosphere with (3 < 1 (see FigO, and 
is much easier to move with j3 > 1. 

In a weak ambient magnetic field {j3 > 1), the cloud ex- 
pands while it moves. The evolution naturally leads to the 
strong thermal insulation of the cloud because it is fast enough 
to drag the magnetic field so as to be "wrapped" by it, and to be 
thermally insulated from the surroundings. For weak magnetic 
fields, the closed field topology yields the best thermal insu- 
lation because the magnetic field perfectly envelops the rising 
cloud. However, for the same reason, it limits the cloud expan- 
sion and motion. 

The observation supports the result of the linear expansion 
of a factor of 3-4 (Ciaravella et al. 2000). In the circumstances 
of the cylindrical clouds, we find an expansion factor F e w 
2.1 for the closed magnetic field simulation, F e = 3.9 for the 
open magnetic field and an intermediate value F e w 2.9 for 
the basic hydrodynamic case. Thus, the open field favours the 
cloud expansion, while the closed field does not. The expansion 
of the cloud nearly stops at the end of the simulations because 
the pressure equilibrium i s reestablished be t ween the cloud and 
the surroundings, and, as lRiley & Crooken (2004) pointed out, 
the magnetic tension is not important as restoring force. 

We find similar results and an expansion factor F e = 3.6 
for the simulation of a cylindrical cloud which carries an initial 
moderate magnetic field component along its axis. A strong in- 
ternal field does not appear realistic in our conditions because it 
would imply an extremely cool cloud all along its evolution, in 
contrast with observations. The effect of the internal magnetic 
field component is to further funnel heat transport and therefore 
reduce thermal exchanges with the surroundings. However, we 
should expect thermal exchanges at the extremes of a more re- 
alistic finite cylindrical cloud. 

We conclude that the open field topology is most conducive 
to both the thermal insulation and a good degree of expan- 
sion, and may therefore be the best match to observations. 
This is true regardless if we consider a purely plasma cloud 
or a cloud accounting also for the ejection of magnetic flux. 
Several CME models already considered an open ambient mag- 
netic field in the outer corona in order to take into account 



a steady solar wind flow ([Manchester et al.ll20 04; Chane et al 
20051 iLugaz et"alll2005l) . Moreover. ICargill & Schmidt! d20Q ' 



considered the model of an emerging flux rope evolving in a 
open magnetic field. 

Our work shows that the cloud mass and shape are also 
important. The expansion factor of an initially spherical cloud 
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increases by ~ 4% more in an ambient open field than with 
no field. Instead, for an elongated cloud the increase is of 
~ 35%. Thus, although the open field induces a larger expan- 
sion regardless the cloud geometry, this effect is more signif- 
icant for an elongated cloud than for a spherical one. This is 
in agreement with many models that consider a filament erup- 



gering of a CME ( 


Forbes & Isenberdll99lt Linker et alJl2003: 


Moore et al. 2001; 


T6r6k& Kliem]2005). 



We finally remark that we assume conditions of ideal MHD 
(except for numerical diffusivity that gives an effective mag- 
netic Reynolds number R m > 50). It will be interesting to 
investigate in the future the role of the magnetic diffusion and 
joule heating, which may be important locally in some areas, 
in which the magnetic field is found to be enhanced during the 
propagation of the cloud. 
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